Ring-shaped spatial pattern of exciton luminescence formed due to the hot carrier 
transport in a locally photoexcited electron-hole bilayer 
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A consistent explanation of tlie formation of a ring-sliaped pattern of exciton luminescence in 
GaAs/AlGaAs double quantum wells is suggested. The pattern consists of two concentric rings 
around the laser excitation spot. It is shown that the luminescence rings appear due to the in-layer 
transport of hot charge carriers at high photoexcitation intensity. Interestingly, one of two causes 
of this transport might involve self-organized criticality (SOC) that would be the first case of the 
SOC observation in semiconductor physics. We test this cause in a many-body numerical model 
by performing extensive molecular dynamics simulations. The results show good agreement with 
experiments. Moreover, the simulations have enabled us to identify the particular kinetic processes 
underlying the formation of each of these two luminescence rings. 

PACS numbers: 71.35.-y, 78.60.-b, 78.67.De 
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INTRODUCTION 

Non-equilibrium collective effects in the exciton and 
exciton-polariton systems in semiconductor heterostruc- 
tures are a subject of intensive studies [l-26|. A par- 
ticular attention has been focused on the beautiful phe- 
nomenon discovered experimentally in the system of in- 
terwell excitons in GaAs/AlGaAs double quantum wells 
(QWs) 0: at sufficiently high excitation intensity, a lo- 
cal photoexcitation of electrons (e) and holes (h) above 
the exciton resonance gives rise to a macroscopic ring- 
shaped pattern of spatial distribution of the exciton lu- 
minescence. The radius of the pattern can be varied in a 
wide range by tuning external parameters such as excita- 
tion intensity or gate voltage. Remarkably, at sub-Kelvin 
lattice temperatures the external ring of the stationary 
pattern exhibits a sharp fragmentation, which could be 
the signature of a non-equilibrium macroscopic quantum 
effect. 

To understand the nature of the ring-shaped pattern, 
one should build a many-body model that captures lo- 
cal generation of electron-hole pairs and their spatial dy- 
namics accompanied by the processes of formation and 
recombination of excitons. If the exciton lifetime is long 
enough, one should also consider the spatial dynamics of 
the excitons. 

The first theoretical explanation which met these re- 
quirements was based (ij, [3, IE Hi on the diffu- 
sive transport model (DTM) applied to the locally photo- 
generated holes and equilibrium electrons which initially 
were uniformly distributed in the quantum-well plane. 
The overlapping region for the hole and electron spatial 
distributions apparently gave rise to the ring of exciton 
luminescence. 

However, this explanation has a lot of evident short- 
comings \27\ . For example, if one adds the same number 
of photogenerated electrons in the model (in principle, 
they must be added to maintain the electroneutrality) , 



then the ring of exciton luminescence can disappear due 
to the exciton formation term, which is simply propor- 
tional to the product of the electron and hole densities. 
(If these densities decrease monotonically from the exci- 
tation spot center, the luminescence intensity would ap- 
parently follow them.) More specific shortcomings can 
be found in Appendix A. 

In this paper a novel consistent explanation of the ring- 
shaped pattern formation is given. The main idea is that 
there appears an essential in-plane electric field in the 
excitation spot region at high enough excitation power. 
This field strongly affects the spatial dynamics of the 
photogenerated electrons and holes. (We do not consider 
any equilibrium carriers at all.) It is shown that there are 
contributions to the electric field from two quite different 
physical processes. Essentially, due to the one of these 
contributions the ring-shaped pattern formation could be 
understood in the paradigm of self-organized criticality 
(SOC) [28|. To test the contribution, we have performed 
extensive molecular-dynamics simulations. They have 
shown that only this contribution is quite enough for 
the detailed qualitative explanation of the ring-shaped 
pattern. (However, this paper does not concern the tran- 
sition to the SOC regime. The parameters for the sim- 
ulations have been chosen to be in this regime from the 
very beginning.) 

This paper is organized as follows. Further in this 
Section we introduce some essential properties of dou- 
ble quantum wells and interwell exciton formation, some 
experimental results we focus on, and the formulation 
of the problem for the research. In Section 2 we sug- 
gest two qualitative explanations ("scenarios") for the 
ring-shaped pattern formation and make the correspond- 
ing estimates. Section 3 describes many-body dynamical 
model and conditions of the molecular-dynamics simu- 
lations performed to investigate the second scenario in 
more details. Section 4 contains the results of the simu- 
lations. Section 5 is conclusion and discussion. Finally, 
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Figure 1: Left: Schematic of a double quantum well (DQW). 
The DQW structure is bordered by highly-doped GaAs lay- 
ers that serve as external electrodes forming a plane capaci- 
tor. In experiments [l^ . [2^ widths of the layers between 
the electrodes are (200 nm){8 nm)(4 nm)(8 nm) (200 nm), 
respectively. Right: (Top) Schematic of (i) the DWQ energy 
profile along the "growth" axis (Z axis) when there is a volt- 
age applied between the external electrodes (it results in a 
linear bias of the profile) and (ii) interwell exciton formation 
(arrows show the path of a photoexcited electron). (Bottom) 
Interwell excitons as co-directed classical dipoles. 
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Figure 2: Experimental curves taken from Ref.Q: lumines- 
cence intensity of interwell excitons vs distance r from the 
excitation spot center at different excitation powers (numbers 
near curves, in nW). Excitation spot radius is about 20 ^m. 
Inset: Dependence of the external ring radius on excitation 
power. 



Section 6 consists of three Appendixes. 

The structure of double QWs used in the experiments 
[^.[isl 24 1 is shown in Fig. [TJleft). The electron band-gap 
energy Egap of the "barrier" (B) layers is larger than Egap 
of the "weh" (W) layers (Fig. fright)) so that GaAs lay- 
ers form two rectangular potential wells with the depth 
Uqw = [Egap (B) ~ Egap {W)) /2 « 0.4 eV . At a mod- 
erate occupation of the wells (i.e., when the number of 
electrons in the GaAs conduction band is not macroscopi- 
cally large; see next Section for details), a voltage applied 
to the external electrodes provides a constant tilt of the 
DQW potential profile (Fig. [Ijriglrt)). This "gate" volt- 
age, Vg, is needed to separate electrons and holes in the 
different wells facilitating the formation of interwell exci- 
tons. The crossover between the interwell and intrawell 
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lace at V„ « 0.3 V ^ 



M^g > 0.3 V, the 



exciton "ground" states takes 
(Since in the experiments 
intrawell excitons are not discussed further.) The sta- 
tionary laser pumping comes along Z axis and is used 
for the formation of a macroscopic number of photogen- 
erated electron- hole pairs. In experiments d, [iS, 24 1 the 
typical laser power was several hundreds of and was 
focused in a spatial spot of few tens of ^m. The pump- 
ing energy was well above the exciton resonances so that 
an electron was photoexcited to a high-energy level of the 
QW near the continuum. Due to the applied gate voltage, 
the electron can then tunnel to the corresponding QW. 
(The effective mass of an electron in GaAs is seven times 
smaller than the one of a heavy hole so the tunneling 
is essentially more probable for electrons than for holes. 
Note that heavy holes have lower energy than the light 
ones in GaAs.) The tunneling might also be facilitated 



by the voltage-induced triangular profile of the barrier. 
Finally, after a spatiotemporal in-layer relaxation of the 
charge carriers, the interwell excitons are formed; they 
live some finite time and then annihilate giving rise to 
the photoluminescence (PL) pattern in the QW (X-Y) 
plane. 

Let us now turn to the experimental results that we 
would like to explain (Fig. [5]) . At small excitation power 
the PL spatial profile practically follows the excitation 
spot (see details in [6]). When excitation power exceeds 
some value a thick ring of luminescence appears near the 
edge of the excitation spot. li\ Ref.fi'l it was already seen 
at excitation power 220 fiW . Hereafter we call this ring 
as "internal ring". Finally, when the excitation power 
exceeds another critical value (Fig. [2]), a thin "external" 
ring of luminescence appears around the excitation spot 
and the inner ring. Everywhere in this paper the words 
"ring-shaped pattern" mean these two concentric rings. 

The formation mechanism of the ring-shaped lumines- 
cence pattern is a subject of the research described below. 
In particular, we pose the following questions, (i) Why 
does the ring-shaped pattern appear only when the laser 
excitation power exceeds some critical values? (ii) What 
are the kinetic processes which underlie the formation of 
the internal and external rings? (iii) Why does the exter- 
nal r^g radius depend strongly on the static gate voltage 
Vg [30|? (iv) Why the luminescence of intrawell excitons 
(i.e.; purely "two-dimensional" excitons) does not exhibit 
the ring-shaped pattern [Qj? 
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TWO SCENARIOS OF THE RING-SHAPED 
LUMINESCENCE PATTERN FORMATION 

In general, we believe that the ring-shaped PL pattern 
appears due to the transport of hot uncoupled electrons 
and holes from the excitation spot at high enough ex- 
citation power. During the spatial spread the carriers 
relax in kinetic energy emitting phonons and can eventu- 
ally form excitons relatively far away from the excitation 
spot ("far away" in comparison with the spot radius). We 
suppose that the hot charge carriers are formed due to (i) 
in-layer electric fields that occur at high pumping power 
in the excitation spot region and (ii) high mobilities of 
the charge carriers in GaAs. We suggest two particu- 
lar mechanisms of the electric field occurrence; they are 
described in details in next subsections. 



First scenario: hot carrier transport induced by the 
external gate voltage 

The first scenario concerns the screening of the gate 
voltage Vg by photogenerated carriers in the excitation 
spot (Fig. [3]) at high excitation power. 



Laser 




Figure 3: Schematic of the electric field distribution along the 
growth axis (Z axis) of DQW structure in the vicinity of laser 
excitation spot at high excitation power. (It is supposed that 
the excitation spot is far away from the DQW edges and ex- 
ternal contacts.) In the spot region (green and red ovals in 
the center) the static electric field induced by external gate 
voltage is curved due to the presence of a macroscopically 
large number of photogenerated charge carriers. The hori- 
zontal projections of the field cause in-layer transport (shown 
by thick arrows) of the carriers from the excitation spot. 



excitation spot is D = E -|- AnP. Here E = EqCz is the 
static uniform electric field generated by Vg and P is po- 
larization of the medium. If n is the density of electron- 
hole pairs in the excitation spot of area S and d is the 
average distance between the pairs in different wells (i.e., 
d is the dipole length), then P = — (nS) (ed) / (Sd) = 
— (ne) and D = (Eq — Anne)ez. In the experiments 
[ 29j the DQW structure was considered as an insulator. 
It means that condition 



Eq ^ At: en, 



(1) 



should be fulfilled in the excitation spot region. In this 
case the gate voltage results in a linear slope of the 
DQW potential energy profile along Z axis on the value 
5U{z)^-eEoz (Fig. fright)). 



However, at typical value Vg ^ IV and ti ^ 10 
[29t one gets cEq w eVg/{2L) - 10^ eV/cm, where 
L = 200 nm is the width of external barrier of the DQW 
structure (Fig. [Deleft)), and Airne'^ ~ 10^ eV/cm. So 
at expected densities n ^ 10^^ ~ 10^^ cm~^ in the exci- 
tation spot at high pumping power, the condition ^ is 
not fulfilled there. It means that Z-axis component of the 
resulting field Ez is essentially dependent on z in the ex- 
citation spot region. More importantly, in this case there 
exists an in-plane component Er of the electric field that 
will push both electrons in one layer and holes in another 
layer away from the excitation spot (Fig. [3]). 



Second scenario: hot carrier transport induced by 
the repulsive in-layer interaction 

If the photogenerated electrons and holes do not leave 
the excitation spot for any reason then the higher the 
pumping power P^x the higher the carriers densities in 
the spot. (Since the excitation is off-resonant, the value 
of the exciton formation time is always larger than the 
time of energy relaxation of the carriers to reach the exci- 
ton transition.) Due to the bilayer geometry, there exists 
a value of excitation power at which the carrier densi- 
ties in the spot reach the values when repulsive in-layer 
Coulomb forces between carriers become stronger than 
the attractive interlayer force. To be more specific, let 
us make an estimate of the interaction strength in the 
excitation spot through the dimensionless interaction pa- 
rameter Tg. At small carrier densities the interaction in 
the spot is dipole-dipole one rather than the Coulomb as 
in the case of a monolayer. In particular, for electron (or 
hole) monolayer 



(eVf)/(ftV("^^')) = r/ae ^ 1/^'. 



(2) 



where ag — / (me^) is the Bohr radius, n is the car- 



rier density in the spot, and r 



-1/2 



is the average 



The electric induction in the bilayer volume at the laser distance between carriers. (Hereafter we omit dielectric 
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constant in the formulae unless making numerical esti- 
mates.) Thus, for monolayer the increase of density n 
leads to the decrease of interaction. However, in the case 
of electron-hole (e-h) bilayer at f > d, where d is the 
interlayer distance, the interaction is dipole-dipole one, 
U = e^dP /f'^ rather than e^/f. This leads to 



r, = C//(nV(mf2)) = (fas) - (d/as) 



(3) 



where n < < ~ lO^^ ^,^-2 ^ ^ iq-Q |29| |. 
It is seen that in this case the interaction increases in 
accord with the densities. Al n — Uc (that corresponds 
to some critical excitation ■power {Pex)^ the character 
of interaction is changing: the repulsive in-layer interac- 
tion becomes dominating and, moreover, one should put 
n = effectively in the estimate ([2]), i.e., the repulsive 
interaction becomes huge. It leads to the appearance 
of in-layer electric fields ejecting the electrons and holes 
from the excitation spot region. Then the e-h densities in 
the spot grow to the critical values again and the ejection 
process recurs. There is a direct correspondence between 
this self-organized ejection, which keeps the critical val- 
ues of carrier densities in the spot , and the avalanches in 
classical sand-pile model of SOC 28l 31|. 

Since the mobilities of charge carriers in GaAs QWs are 
very high (up to 10'' cm^/ (Vs) for electrons at sub-Kelvin 
temperatures [32-35|), the initial kinetic energies of some 
part of ejected particles can exceed the threshold of op- 
tical phonon emission. Their relaxation in energy is then 
very fast so these carriers will likely form excitons not far 
away from the excitation spot. In turn, the carriers with 
kinetic energy less than optical phonon energy fvjjopt go 
further. These carriers relax relatively slow by emitting 
acoustic phonons. In GaAs fvjJopt ~ 37 meV so the veloc- 
ities of the carriers contributing to long-distance trans- 



port are less than Vn 



10 cm/s. (Some more details 



about carrier energy relaxation due to phonon emission 
can be found in Appendix B.) 

In addition, there is an essential difference in mobili- 
ties, and in effective masses, for electrons and heavy holes 
in GaAs 32h34|. This can lead to some difference in 
in-layer "stream" velocities of electrons and holes. The 
latter, in turn, would result in a suppression of the at- 
tractive interlayer Coulomb force (see Appendix C for 
details) until the velocities become small enough due to 
the acoustic phonon emission. In these conditions the 
interlayer exciton formation is also suppressed in some 
region of distances from the excitation spot. (Recall, for 
example, that the classical scattering cross-section for the 
Coulomb potential, the Rutherford cross-section, is pro- 
portional to T^""*, where V is the initial relative velocity 
at infinite distance.) The suppression leads naturally to 
the formation of a luminescence ring that would define a 
circumference of the total luminescence pattern. 

Thus, the internal ring of luminescence (Fig. [5]) might 
appear due to the electrons and holes that had emitted 



optical phonon(s) and then quickly formed excitons. In 
turn, the external luminescence ring can appear due to 
the carriers that were under the optical-phonon emission 
threshold so they needed larger time (and corresponding 
distance) to relax emitting acoustic phonons. 

In general, we suggest that at high photoexcitation 
power there exists an in-plane electric field which 
consists of two contributions, "gate voltage-induced" and 
"in-layer interaction- induced". Due to high mobilities of 
charge carriers in GaAs even a moderate value of Er re- 
sults in high initial velocities of the carriers directed to- 
wards the outside of the excitation spot. Large relative 
velocities of the ejected electrons and holes lead to the 
suppression of the interlayer Coulomb attraction between 
them. The latter results in the suppression of exciton for- 
mation in some domain of distances from the excitation 
spot. In turn, the formation dynamics of ring-shaped 
luminescence pattern can be divided in three stages: (i) 
radial acceleration of carriers in the excitation spot region 
due to the in-plane component of the static electric field 
that appears at relatively high carrier density in the ex- 
citation spot and due to the in-layer Coulomb repulsion 
at high pumping power; (ii) slowing down of unbound 
carriers due to emission of optical and acoustic phonons 
and due to the ambipolar electric field ("Coulomb drag") 
(iii) the regime of strong interlayer Coulomb correlations 
- formation and optical recombination of interlayer exci- 
tons. 

In what follows we focus on the second scenario and 
test it by molecular dynamics (MD) simulations. 



MOLECULAR DYNAMICS SIMULATIONS: 
NUMERICAL MODEL 



To describe the spatial dynamics of N hot electrons 
and holes, we use classical equations of motion 
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combined with the conditions of exciton formation and 
optical phonon emission (see below). Here vectors r* 
and are in-plane positions of i-th electron and j-th 
hole (1 < i,j < N), fn*f,!^h) is the electron (hole) effective 
mass, e is the electron charge, and d is the interlayer 
distance (see Fig. 2]). 

In addition to the inertia terms, the left-hand side of 
Eqs.(|3][S]) contains phenomenological momentum damp- 
ing terms due to the interaction with acoustic phonons 
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Figure 4: Qualitative schematic of optically-excited electron- 
hole bilayer. Both stationary laser pumping in the center and 
spatially-distributed luminescence are perpendicular to the 
layers. The charge separation between the layers postulated 
in the numerical model is due to the external gate voltage 

Vgate • 



with constants "fe{h) = e//ie(/i), where fj.e{h) is electron 
(hole) mobility. The dimensionless equations read 
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(6) 
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with constants ci = m*fj,e/ (m'^iifi) ^^nd C2 = ml/m'^. 
Hereafter, we normalized time by = •\/e'7i*/ie/e, where 
e is dielectric constant of the layers, and all distances 
by = \/rne/Zg. To estimate the parameters, we used 
well-known experimental values for high-quality undoped 
GaAs/AlGaAs QWs. In particular, taking typical /ig ~ 
10"^ cm2/(Vs) [si, [sil for temperatures T < 1 K, m* sa 
0.067me, m'^ « 0.5me (we is bare electron mass), e = 
12.8 and Hh ^ O.l/Xe 0, one gets ^ 10"^ s, ^ 10^^ 
cm, ci ^ 1 and C2 ^ 0.1. 

The optical phonon emission was modeled in the fol- 
lowing way: if the kinetic energy of a carrier exceeded 
the energy of optical phonon, the latter was subtracted 
from the first and the new direction of the carrier velocity 
became random. 

Simulation of the Laser Pumping. Stationary op- 
tical pumping of carriers was simulated by generating 
them in random positions inside the excitation spot of 
radius ro with some generation rate which was modeled 
in two different ways. In fact, during a MD simulation 
the time is changed by discrete steps, with the elemen- 
tary time step At. According to the first way [isl [36|. 
the generation rate p was defined as a probability per 
At to create one e-h pair in the excitation spot, so that 
pAt < 1. We refer this case as single generation regime 
( SGR ). Alternatively, one can consider the formation of 



several e-h pairs during At. Then the carrier generation 
rate (CGR) is defined as a number of e-h pairs generated 
in the excitation spot during the time step At. We call 
it as multiple generation regime (MGR). Note that the 
results of MD simulations differ essentially in the single 
and multiple regimes. Indeed, it is intuitively clear that 
SGR likely corresponds to weak pumping whereas MGR 
describes high-power excitation. 

The initial velocities of carriers in the excitation spot 
were also chosen randomly within intervals | r* | < and 
|i";J < 'ywo, where we took = 0.5 for all simulations. 

During the spatial dynamics of the carriers, the exci- 
ton formation happened if an electron and a hole were 
close enough to each other, |re — r/i| < a, where a{d) 
is a phenomenological in-layer exciton radius, and their 
relative velocity was smaller than some critical value, 
|re — rft,| <Vc [25| (see also Appendix C). Note that the 
dependence of exciton formation rate on the e-h rela- 
tive velocity is one of the most crucial ingredient for the 
ring-shaped pattern formation: assuming the absence of 
such dependence one always gets a spatially-monotonic 
decrease of the luminescence outside tlie excitation spot 

To simplify the simulations, we did not consider the 
exciton dynamics. It means that as soon as an electron 
and a hole had formed an exciton, their dynamics was no 
longer considered and the position of the formation event 
was recorded as a position of photon emission. Qualita- 
tively, this corresponds to zero exciton lifetime. 

One should note that since both the internal and ex- 
ternal ring radii are temperature- independent [2^ , it 
is not advisable to include temperature (i.e., to add a 
stochastic force in Eqs.(|3]lS])) in the consideration. In 
turn, since the low-temperature fragmentation of the ex- 
ternal ring apparently depends on the exciton dynam- 
ics, we do not expect to observe the fragmentation in the 
simulation results. 

Finally, due to the inevitable restrictions in computa- 
tional power it was only possible to simulate the dynam- 
ics o{ N < 10'' interacting particles. For this reason we 
had to modify the values of a, vq, Vc, d etc in compari- 
son with realistic values to facilitate the exciton forma- 
tion. However, it was clearly seen that the closer values 
of those model and real parameters, the better the corre- 
spondence between the MGR simulation results and the 
experimental ones. Note that the in-plane motion of the 
carriers was not restricted by any spatial boundaries. 



MOLECULAR DYNAMICS SIMULATIONS: 
RESULTS 



Some preliminary results for the SGR were published 
in Ref.pS]. In particular, quasi-lD simulations and the 
crucial dependence of the ring pattern formation on the 
critical relative velocity Vc in the exciton formation con- 
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dition were discussed there. In what follows, some new 
essential results are described. 



same behaviour, i.e., independence of the ring position 
on p, was observed in quasi-lD case (25| . 



Single generation regime (SGR) 

The results of MD simulations of Eqs.dMT]) in the SGR 
for several sets of parameters indicate that there are two 
qualitatively different pictures. In general, the in-layer 
distribution ni^^ir) of stationary luminescence exhibits 
a ring-shaped pattern around the excitation spot. How- 
ever, the pattern always contains only one ring. More 
importantly, the ring can originate by two qualitatively 
different ways. 

According to the first way, the in-layer distributions 
of electrons and holes are separated from each other and 
the ring occurs in the region of their overlapping (Fig. 
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Figure 5: Stationary in-layer distributions of electrons (rie), 
holes (rifi) and luminescence (nium('")) averaged over time in- 
terval (200,830). Inset: dependence of the total number A'^ of 
electrons and holes on time t. In a stationary state the car- 
rier generation rate is balanced by the luminescence rate so 
N{t) exhibits saturation. Parameters of the simulation: exci- 
tation spot radius ro = 4, p = 10, velocity of optical-phonon 
emission threshold Vopt ~ 50, maximal electron initial veloc- 
ity vo = 50, critical relative velocity Vc = 10, critical relative 
distance a = 0.2, first and second coefficients in the equations 
of motion for holes ci = 1 and C2 = 0.25, = 0.01, time step 
At = 0.0005. 



The dependence of the ring position on generation rate 
P (Fig. Oj which mimics the excitation power, shows that 
though the luminescence ring intensity increases with p 
its position is virtually independent on p (Top inset in 
Fig. [6]). This behaviour differs from that observed ex- 
perimentally (Fig. (2), where the radius of external lu- 
minescence ring grows near linearly with the increase of 
photoexcitation intensity and the growth of the internal 
ring radius is also quite noticeable. Interestingly, the 
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Figure 6: Dependence of stationary in-layer luminescence dis- 
tribution niuni(r) on generation rate p. Top inset: dependence 
of the luminescence ring radius on p. Bottom inset: stationary 
total number Nstat of electrons and holes vs p. Parameters of 
the simulations are the same as in Fig. [5] 



However, there exists an another way of the ring pat- 
tern formation. It was observable at other sets of param- 
eters, in particular, when maximal initial velocities of the 
carriers, the critical relative velocity Vc and the distance 
a were relatively small. Note that the CGR p was taken 
in the same range of values as previously. 

According the second way, the in-layer distributions of 
electrons and holes practically coincide with each other 
(Fig. [T]). The ring occurs at the outer side of the dis- 
tributions, where the differences in the carrier densities 
and velocities are small enough to allow the exciton for- 
mation. The dependence of the ring radius on the CGR 
p (Fig. [5]) has shown that, in contrast to the previous 
case (Figs. [5] and [5]), the radius increases linearly with p. 

Summarizing the simulation results for the SGR, one 
can conclude that the first way of the ring pattern for- 
mation, when electron and hole in-plane distributions are 
separated and the ring is formed at their overlapping, 
does not correspond qualitatively to the experimental re- 
sults d, [H, Q. The second way could mimic the ex- 
perimental situation when the excitation power was such 
that only the internal ring of luminescence was observ- 
able. However, the position of the ring in the simula- 
tions depends on the pumping rate stronger (linearly) 
than that of the internal ring (cp. Fig. [5] and Fig. [5]) in 
the experiments 0, [l3, 24|. In turn, nearly- linear depen- 
dence of the ring radius on the excitation power is typical 
for the external ring of luminescence (Fig. [2]), but then 
the simulations omit the internal ring. 

This misfit can indicate that though the parameters 
of the simulations in the second case are closer to re- 
alistic ones, one needs more realistic carrier generation 



7 




Figure 7: Stationary in-layer distributions of electrons (n^), 
holes (rih) and luminescence (nium(»")) averaged over time in- 
terval (500,10000). Top inset: dependence of the total num- 
ber A*' of electrons and holes on time t exhibits a station- 
ary state. Bottom inset: dependence of electron (14) and 
hole (Vh) stream velocities vs distance r from the excitation 
spot center. Parameters of the simulation: ro = 3, p — 10, 
Vopt = 10, uo = 1, 14 = 1, a = 0.05, ci = 1.1, C2 = 0.3, 
= 0.9 (at smaller cf the results are essentially the same), 
time step At = 0.005. 




Figure 8: Dependence of stationary in-layer luminescence dis- 
tribution niuin('') on generation rate p. Top inset: dependence 
of the luminescence ring radius on p. Bottom inset: stationary 
total number Nstat of electrons and holes vs p. Parameters of 
the simulations are the same as in Fig. [T] 



algorithm which would enable us to model high excita- 
tion powers and, at the same time, is independent on the 
MD time step At. 



Multiple generation regime (MGR) 

In the MGR the simulations result in much better cor- 
respondence with the experimental plots (Fig. [5]) and, 
simultaneously, one can trace some correspondence with 
the SGR results ("the second way"). 
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Figure 9: Top: Stationary in-layer distributions of electrons 
(ne), holes (uh) and luminescence (niuni('')) averaged over 
time interval (100,400). Inset on top: dependence of the total 
immber A'^ of electrons and holes on time t exhibits stationary 
regime. Bottom: 3D luminescence pattern. Parameters of 
the simulation: ro — 3, pumping rate (number of particles 
generated during time step At in the spot) ANapot = 3, Vopt = 
10, «o = 1, 14 = 1, a = 0.05, ci = 1.1, C2 = 0.3, (f = 0.9, 
At = 0.01. 

First, the MGR results show two concentric rings in 
the in-layer luminescence pattern: the internal ring has 
small radius and high intensity and the external ring has 
relatively large radius and is weaker in intensity (Fig. 

The in-layer distributions of electrons and holes are 
similar to those of the above SGR "second way" (Fig. 
[7]), i.e., the distributions practically coincide with each 
other. 
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Due to the increase of effective pumping power in the 
MGR case the total number of carriers in the stationary 
state (inset in Fig. ^ is more than one order of magni- 
tude larger than that in the SGR case. Note that the sta- 
tionary number of carriers in the excitation spot is very 
close to the product of instant CGR {ANgpot/'^t — 300) 
and the unit time (=1). It means that most of the 
carriers escape from the spot very fast. This can be 
understood by making an estimate of the critical num- 
ber of carriers in the spot required for dominant in-layer 
repulsion (see Section 2): taking ric ^ one gets 
Nc = ?t-cTq 7r(ro/(i)^ ~ 30, which is much smaller 
than the above number of carriers formed in the spot 
during the unit time. (For the estimate we have used 
the parameters in Fig. [9] caption.) Thus, the system of 
charge carriers is in the critical state. 

The mechanism of the ring pattern formation in these 
conditions is as follows. The internal ring is formed due 
to the carriers which have emitted optical phonons. It 
can be seen in Fig. [10] by the sharp changes of carrier 
velocities within the excitation spot. Note that the ve- 
locities essentially exceed the maximal initial velocities 
there. (In addition, the optical phonon emissions were 
seen during the simulations by changes of instant maxi- 
mal velocities of the carriers that were up to the thresh- 
old velocity of optical phonon emission.) It means that 
the in-layer repulsive Coulomb interaction is dominant 
in the excitation spot. The repulsive forces accelerate 
the carriers in such a way that the major part of them 
emits optical phonons and then quickly forms excitons. 
The carriers with velocities beneath the threshold veloc- 
ity of optical phonon emission go further emitting acous- 
tic phonons and, eventually, form excitons relatively far 
away from the excitation spot. This description is in ex- 
cellent agreement with all sets of simulations in the MGR 
so it justifies the second scenario suggested in Section 2. 
Some details of the Coulomb repulsion accompanied by 
phonon emission are given in Appendix B. 

To prove the crucial role of Coulomb interactions, we 
have performed the simulation where the interactions are 
switched off even though the exciton formation condition 
holds. The results are shown in Fig. [TlJ One can see 
that the saturation of N{t) (i.e., the stationary state) is 
absent during the time which essentially exceeds previ- 
ous simulation times for the MGR. The in-layer distribu- 
tions of electrons and holes resemble those for the "first 
way" SGR, i.e., the distributions are spatially separated. 
However, the ring does not form and the luminescence 
decreases monotonically from the center. 

We now discuss the dependence of the ring-shaped pat- 
tern on the pumping power and the critical relative ve- 
locity Vc from the exciton formation condition. 

The dependence of in-plane positions of the rings on 
the pumping rate is shown in Fig. [12] The external ring 
radius increases nearly linearly whereas the internal ring 
radius grows more slowly. This behaviour exhibits a good 
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Figure 10: Dependence of electron (14) and hole (Vh) stream 
velocities vs distance r from the excitation spot center. One 
can see that the carrier velocities within the excitation spot 
region essentially exceed the maximal initial velocity (= 1). 
Note that 14 and 14 are the velocities of that part of the 
carriers which are below threshold of optical phonon emission. 
Inset: Electron (/e(l^)) and hole {fh{V)) distributions over 
single-particle velocities. Parameters of the simulation are 
the same as in Fig. [O] 



agreement with the experimental curves (Fig. (To 
determine the dependences more accurately one needs to 
collect larger statistics that is a very time-consuming pro- 
cedure due to the large total numbers (^ 10*) of particles 
in the MGR.) 

The dependence of the ring-shaped luminescence pat- 
tern on the critical relative velocity 14 shows (Fig. [T5)) 
that the smaller Vc the larger the external ring radius 
(top inset in Fig. [T5|) and the ring intensity. At the same 
time the stationary number of carriers grows up to the 
limiting values of order of lO** (bottom inset in Fig. [T5|) 
for the available computational power. 

Finally, we have performed simulations with two iden- 
tical but spatially separated excitation spots to com- 
pare our results (see Fig. I14p with experimental pictures 
[ll,[33|. One can see that when the spots are placed close 
enough, the external rings open towards each other form- 
ing a figure similar to one of the Cassini ovals. This be- 
haviour also corresponds to the experiments even though 
the simulations do not include the exciton dynamics. 

Summarizing the simulation results for the MGR, we 
conclude that they show reasonable correspondence with 

2ia 




Therefore the theoreti- 
leRll m of the ring pattern 



the experiments 
cal explanations 
formation based on the diffusion-induced spatial separa- 
tion of in-layer distributions of electrons and holes (see 
Appendix A) must be revised. 
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Figure 11: In-layer distributions of electrons (rie), holes 
{uh) and luminescence (niumir)) averaged over time interval 
(1500,2700) in the absence of any direct Coulomb interac- 
tions. Inset: dependence of the total number A*' of electrons 
and holes on time t. It shows that the dynamics is still non- 
stationary. Parameters of the simulation are the same as in 
Fig. [3 



3000 




Figure 12: Dependence of stationary in-layer luminescence 
distribution nium{r) on excitation power with the number 
ANspot of particles generated per time step in the spot shown 
by numbers (1, 3, 6) near the corresponding curves. Inset: 
dependence of the external ring radius on excitation power. 
Other parameters of the simulations are the same as in Fig. 

m 




Figure 13: Dependence of stationary in-layer luminescence 
distribution niuni('') on critical relative velocity Vc from the 
exciton formation condition. Top inset: dependence of the 
external ring radius on Vc- Bottom inset: stationary total 
number Nstat of carriers vs Vc- The pumping rate ANspot = 
3. Other parameters of the simulations are the same as in 
Fig. [3 




Figure 14: Left: Luminescence patterns for the case of two 
excitation spots at two distances, 85 (top) and 70 (bottom), 
between the spots. Right: The corresponding dependences of 
total numbers Ni^2 of carriers excited by the spots on time 
t show that the system arrives at the stationary state. Pa- 
rameters of the simulations are the same as in Fig. [9] except 
critical relative velocity, here Vc = 0.5. 



CONCLUSION AND DISCUSSION 

It has been shown that the stationary ring-shaped lu- 
minescence pattern forms due to the hot carrier transport 
(HCT) caused by the in-plane electric fields which appear 
at high enough excitation power in the excitation spot re- 
gion. The HCT is essentially non-diffusive. In particular, 
the internal luminescence ring appears due to electrons 



and holes emitting optical phonons, whereas the external 
ring forms due to the relaxation of the carriers that are 
initially below the threshold of optical phonon emission. 
To form excitons, these carriers relax emitting acoustic 
phonons and, in addition, due to the interlayer Coulomb 
drag. 

The ring-shaped pattern formation is particularly in- 
teresting as a possible bright signature of self-organized 
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criticality |28l |31 |. Though "the second scenario" nat- 
urally involves the SOC regime, the MGR simulations 
reported have been performed in the critical state. So 
the transition to the SOC regime as well as its statistical 
properties (e.g., 1//- noise) in this system are still open 
questions. 

The author thanks L. P. Paraskevova and L. I. Kon- 
drashova for the encouragement, and Yu. M. Kagan and 
F. V. Kusmartsev for helpful discussions. 



APPENDIXES 

Appendix A: Diffusive model of charge carrier 
transport 

The diffusive transport model [13] used to explain the 
experiments |^ was based on two reaction-diffusion equa- 
tions 

rie = DeV^rig - wneUh + Je (r) , (8) 
rih = DhV'^nh - wUeUh -f Jh (r) , (9) 

where rig and Uh are electron and hole two-dimensional 
(2D) densities, w is electron-hole binding rate to form 
an exciton. The source term Jh (r) = PexS (r) for pho- 
toexcited holes is focused in the local excitation spot. 
The density of photoexcited electrons is supposed to be 
negligible in comparison with equilibrium electron den- 
sity rioo in the absence of laser excitation. When rioo is 
spatially disturbed due to the presence of holes, there 
appears the electron current Je (r) = I — a - Ue (r) , which 
is spread in the quantum-well plane. Here, / and a ■ Ue 
are the currents in and out of the system, respectively, 
such that rioo = I /a- (Note that a here is not the critical 
relative distance used in the simulations but an indepen- 
dent parameter.) Implying the stationary regime and 
the symmetry over the polar angle, and neglecting the 
exciton diffusion [l^, one gets the exciton PL intensity 
JpL(r) (X ne.(r)nh{r). 

Authors of Ref . [13| have assumed that a luminescence 
ring with radius R appears at the overlap of the electron 
and hole densities (see Fig. [T5)) so that rih » Ue at 
r < R and n/j << rig at r > i? with 71^ (r oo) = riao 
(Fig. [15]). Neglecting the exciton formation term wneUh 
far from the boundary r = R, one gets for holes 
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Figure 15: The density distributions (the scales are in a.u.) 
obtained in Ref. [13] with the use of the diffusive transport 
model described in text. 



with boundary conditions Tig (r = i?) = 0, (r — > oo) ~ 

Denoting A = (ng — tIqo) /?^oo and x — r/X, where 
A = ^ Del a — \J Detiool I is a characteristic length, we 
arrive at 



1 dA 

X dx 



A = 0. 



(12) 



The last equation is the modified Bessel equation of zero 
order. The solution is 



Aix)^A-Ko (x) , 



(13) 



where A is some constant and Kq {x) is modified Han- 
kel (or MacDonald) function of zero order, such that 
Ko (x ^ 0) « In i, Ko (x^co) ^ 0. Thus, 



Tie (r) /noo = l + A-Ko (r/A) . 



(14) 



If A ^ i?, as supposed in [l3j, then at i? < r << A the 
electron density (r) /rioo ~ 1 + Aln(A/r). Using the 
boundary condition Ue (r = R) =0 one finds coefficient 
A. Finally, 

Tie (i? < r << A) w [1 - In (A/r) / In (A/i?)] . (15) 

Then it has been assumed that at the boundary between 
the electron and hole densities the total current is zero, 



Using boundary condition Uh {r = R) ^ it results in 



rih (r < R) 
For electrons one gets 



2TTDh V r 



(I/De) + (a/De) ■ Ue (r) 



(11) 



One. 



From here it follows that 

Denoo/\n{X/R)=Pex/2TT, 



(16) 



(17) 
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so the ring radius can be expressed as 

R = Aexp {-2t: DeUoc/Pex) ■ 



(18) 



If we set Dfi = 0, then the ring radius R must be equal 
to zero since, according to the model [l3|, the diffusion 
of holes is the only reason why they move out of the 
excitation spot. However, Eq. p8p does not depend on 
Dh so the ring radius is not zero at Dh = 0, i.e., the ring 
exists even if all holes are left into the excitation spot. 
This clearly unphysical result is not a consequence of the 
limiting case R, but rather comes from wrong initial 
assumptions. (In Ref.flJ] the result, Eq. ([T8|) . has been 
generalized, but even then R does not dep end on £*/,..) 

We note that in the original paper [13| the erroneous 
formula for R has been given |27| : 



Rorig = A exp (- 



2TTDenool {DhPex)) ■ 



(19) 



The presence of Dh in the exponent denominator could 
be deceiving since at first sight (i.e., without dimen- 
sionality check: [D^ — [Dh\ = cw? js^ [n-oo] — cm~^, 
[Pex] — s^^) it looks reasonable. 

In addition, both the DTM |13| and its modifications 
0,116-18] could not explain in principle why the exter- 
nal luminescence ring appears only when the excitation 
power exceeds some critical value. 

Nevertheless, the drift-diffusion regime can be appli- 
cable for slow charge carriers near the luminescence ring 
at r ~ i?. (As before, we do not consider equilibrium 
carriers and are only focused on photogenerated ones.) 
In particular, the continuity equations in this regime are 
given by 



ne{h) + divie(h) = 5e(/i) ~ r, 
fix +divix = r - ux/tx- 



(20) 
(21) 



Here rig, rih and ig = — ng/ieE — ZJeVrie, ih = 
nh^h^—DhVrih are 2D densities and particle flux den- 
sities of uncoupled electrons in plane z = d/2 and holes 
in plane z = —d/2, He{h) is electron (hole) mobility. The 
particle flux density for excitons ix ~ —Dx^nx, where 
nx is the interlayer exciton density. The contribution 
from the dipole-dipole interaction between the excitons 
is omitted in ix since it appears as an above-linear cor- 
rection on nx ■ The carrier generation rates 5e(/i) (r, t) are 
some given functions. The exciton formation rate can be 
written as (hereafter inessential constant prefactors are 
dropped) 

r (r, t) = Jw (|vi - V2I) fe (r, vi, t) fh (r, V2, i) d^^id^^2. 

(22) 

where (r, v, t) is the electron (hole) distribution 
function, so that ne(;j) (r, t) = //^(^) {r, v, t) d^v, and 
w {v) is the specific exciton formation rate. The exci- 
ton lifetime tx is supposed to be density- independent. 



Finally, the Poisson equation for the electric field reads 
(time dependence is dropped; e is dielectric constant) 

div (eE(r, z)) = 47re[(nft(r) + nx{r)) 5{z + d/2)- (23) 
-{n,{r)+nx{r))5{z-d/2)]. 

It includes the contribution of the interlayer exciton 
dipole fields and keeps the electroneutrality for the free 
carrier system when the exciton formation is suppressed 
(nx(r) = 0). 

At r ^ i? one can put w {v) w ifmax, then F (r) 
WmaxKe {r) Uh (r) and the Eqs.^,^,^ with g^(^h) = 
become a closed system. 

We note that the ambipolar electric field E might play 
an important role in the formation of a sharp intensity 
profile of the external luminescence ring. In this regard, 
it is useful to note that the FWHM of the external ring 
intensity is almost independent on the ring radius R at 
high excitation powers (see Fig. [2]). 



Appendix B 

If \pn > hojopt at n > d~^ the in-layer Coulomb 
repulsion in the excitation spot might be "exhausted" 
at small distances: the potential energy of the carriers 
transforms into kinetic one which, in turn, is spent for 
the fast optical phonon emission so that the carriers do 
not go far from the excitation spot. 

To estimate whether the values of carrier densities in 
the excitation spot are sufficient for such process, let us 
consider two electrons resting at distance vq from each 
other at moment t = 0. We neglect the energy dissipation 
due to acoustic phonon emission first, to facilitate the 
effect described above. Then the equation of motion is 



mr 



2 I 2 
/r 



with m = 'Tt*/2, r = |ri — r2|. (Dielectric constant e is 
introduced by — > e^/e in the final expression.) The 
solution is expressed through the inverse function, 

t/to — \/ — X + i In + \/ — , 



where x = r/ro, io = \J'mr%/ (26^). At t > (3 4)to 
with a good precision the velocity v k, Voo — To/to, i. e., 
the electrons move near uniformly at large times. The 
condition of optical phonon emission reads 

mv'^/2 = (? /tq - (? /r> hiOopt-, 

from where it follows 

~ " \ - hioopt/ {e^ /tq) 
Substituting tq ~ in the expression for Tc, one has 

n ^ (ehuopt/e^ + i/rcY 
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The value n* = {ehujopt/e^ 



4 • 10" cni-2 (e « 12.8, 



fujJopt ~ 37 meV for GaAs) at rc = oo is the smallest 
density at which the process of optical phonon emission 
is dominant. Note that at n > n* the effect is extremely 
pronounced, e.g., at n = l.OOSn* (extra 0.3% to n*) one 
gets Tc ~ 1 /.tm. At n = n*, we obtain tq O.las for 
GaAs. 

To find a qualitative dependence of carrier flux velocity 
on in-plane coordinates when the carrier kinetic energies 
are below the threshold of optical phonon emission, we 
consider the previous model adding the dissipation due to 
acoustic phonons. Then for two electrons we have (unit 
vectors for the Coulomb force are dropped) 



m*ri = e^l |ri - r2 
m*r2 = -e^/ |r 



1^1, 



1-2 



where 7 = e//Ze is dissipation coefhcient and is 
electron mobility. Substituting R = (ri-|-r2)/2 and 
r = ri — r2, in the center-of-mass frame (R — 0), the 
equation of motion reads 

2/2 • 

mr ~ e /r — -^If- 

The dissipation of energy E{r) — mv'^ (r) /2 + e^/r is 
given by dE/dr — ~^v{r) /2 (since dE/dt = — 7w^/2) 
that leads to equation 



dv 1 
mv(r) — = —~ --/V (r) 
dr 2 



In the dimensionless form (r = vqx, v = Woou), we have 



u{x) 



du 
dx 



1 

2^' 



(24) 



where parameter a — y/ 7^rg/ (8me^), x > 1 and u(x = 
1) = 0. (To avoid a confusion, note again that this a 
has its own meaning.) Using relation 7 = e//ie, one 



can rewrite the parameter as a = ^(ro/^e)^^^, where 
= ^m*ii1 is characteristic length scale. Taking for 
GaAs ^ 1 /^TO (see the description of the numerical 
model in the main text) and tq ~ Uc ~ d ~ 10~^ cm, 
one gets a ~ lO^'^. The equation (|24l) can be reduced to 
the Abel equation of the second kind. Its exact analytical 
solution is unknown. However, a qualitative behaviour of 
u{x) can be found from asymptotic solutions at x — >■ 1 
and a; — > 00 (see also numerical solutions in Fig. I16|) . It 
is natural to assume that velocity u{x = 00) = 0. So 
du/dx can be dropped in Eq. (j24p at a; — )■ 00 that gives 
Uqo{x) w l/(2aa;^). Let us now rewrite Eq. as 



du 
dx 



1 



2x'^u 



At X — )• 1, due to the initial condition u — so now the 
last term in the right-hand side can be dropped. This 



0.8 
0.6 

u(x) 

0.4 
0.2 









aToOOl 






^\0.005 


\ \ \ 
\ \ \ 




U.U1 




\o.02 




\ \o.05 






1 \o.i \ 






Vq-sV Y_ 







1 20 40 X 60 80 100 




200 400 ^ 600 



800 1000 



Figure 16: Numerical solutions of Eq. (I24|l with u{x) = v/v^o 
and X — r/ro for different values of parameter a. Top: Linear 
scale. Bottom: Log scale. 



gives Ui{x) « ^yl — l/x at x « 1. The dependence of 
u{x) on a is as follows: the larger a the narrower the 
width of the peak and the smaller its height (see Fig. 
[T6| . Note that in the laboratory frame v^.h = v/2 and 
reM = r/2. 

It is also interesting to note that the MGR simulations 
indicate that the external luminescence ring appears only 
if the carrier stream velocities 14(^)(r) as functions of 
distance r from the excitation spot (see Fig. [TO)) have 
a non-monotonic dependence similar to that in Fig. 1161 
i.e., at small r the velocities grow, reaching the maximal 
values at some finite r (often within the excitation spot 
region), and then they decrease to zero at large r. 
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Appendix C: Suppression of exciton formation in 
the bilayer at high e-h relative velocity 

Here, we illustrate how the critical relative velocity Vc 
appears in principle in the exciton formation condition. 
We can distinguish two mechanisms referred further as 
"geometric" and "kinetic" ones, which lead to the exis- 
tence of the critical e-h relative velocity above which the 
interlayer exciton formation is strongly suppressed. 

1. "Geometric" mechanism. Due to the bilayer 
geometry, at sufficiently large relative e-h velocity the 
interlayer Coulomb attraction, which results in an exci- 
ton formation, is suppressed. To show this, let us con- 
sider the Fourie transform Uq = J cPv exp (iqr) U (r) of 

the pair interaction potential U (r) = — e^/ y/r'^ + d? be- 
tween an electron from one layer and a hole from another 
(r is in-plane relative distance). One gets 



lar to the thread are given by 



Uq = -27re^ 



Jo [qr) rdr 



2^ 



cxp(-(7d), (25) 



the "screened Coulomb potential" in the momentum 
space. (The effects of charge screening studied in the 
random phase approximation ^37i] lead to a change in the 
preexponential factor, which is not important in this con- 
sideration.) Thus, if the electron-hole relative velocity 
V = hq/m > Vc = h/ (md), the interaction between the 
carriers decreases exponentially with the increase of V. 
This means that one can neglect the interaction as well 
as the exciton formation atV>Vc- At d sa 10~^ cm (29| 
and reduced e-h mass m w 0.06me in GaAs, we obtain 
Vc ~ 3 • 10^ cm/s, which is of the same order of mag- 
nitude as the threshold velocity Vmax of optical phonon 
emission in GaAs {vmax/Vc ~ 1.5). 

2. "Kinetic" mechanism. The second mechanism 
is based on the fact that to form an exciton the unbound 
electron-hole pair must emit acoustic phonon. (Here we 
suppose that the carrier velocities are below the threshold 
of optical phonon emission.) 

To illustrate this, let us consider a model system: an in- 
finite train of electrons separated from each other by dis- 
tance L uniformly moves with velocity V along a thread 
and an immovable hole is located at distance d from the 
thread. The interaction potential between the electron 
train and the hole as a function of time reads 



Uit) 



+00 2 



E 



Pk 



where pk = \J (kL + Vt)'^ + d"^. Although the sum di- 
verges as 1/ |fc|, the relative value of the potential (5f/(t) — 
U{t) — t/(0) is convergent. The components of the corre- 
sponding force acting on the hole along and perpendicu- 



{kL + Vt) 



k—~-oo 



pI 



+ 00 



k— — oc 



Both the potential and the force are periodic functions 
of time with period T = L/V, e.g., SU{t + T) = 6U{t). 

Ill the 2D case, when there exists a relative fiow (with 
velocity V) of electrons in one layer and holes in another, 
the interlayer interaction potential between the electron 
ffow and a given hole oscillates with frequency uj ^ V ypn^ 
where n is 2D density of electrons in the ffow. If this fre- 
quency is higher than r~^, where t = Tcivaije-ac^T'h-ac) 
is the minimal carrier-acoustic phonon scattering time, 
the exciton formation in real space is not possible. Thus, 
one can write the condition of exciton formation as 



V < Vr 



(26) 



If the carrier densities are essentially different, one should 
take n = max(ne,n?i) in (|26p . At n ^ 10^" cm~^ and 
T ^ 10~^ s the critical relative velocity is Vc ^ 10"* cm/s. 

In fact, one can get the estimate (pS)) in a more simple 
way. We suppose for definiteness that locally Ue > fih- 
Then an e-h pair with relative velocity V can form an 
exciton if 



Vt < rc 



-1/2 



that is the electron or the hole should have time to emit 
an acoustic phonon before the next electron would come 
to the hole. The condition (|26|) follows directly from the 
last formula. 
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